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Abstract 

We study the formation of cosmic strings by confining a stochastic magnetic field into flux tubes 
in a numerical simulation. We use overdamped evolution in a potential that is minimized when 
the flux through each face in the simulation lattice is a multiple of the fundamental flux quantum. 
When the typical number of flux quanta through a correlation-length-sized region is initially about 
1, we find a string network similar to that generated by the Kibble-Zurek mechanism. With larger 
initial flux, the loop distribution and the Brownian shape of the infinite strings remain unchanged, 
but the fraction of length in infinite strings is increased. A 2D slice of the network exhibits 
bundles of strings pointing in the same direction, as in earlier 2D simulations. We find, however, 
that strings belonging to the same bundle do not stay together in 3D for much longer than the 
correlation length. As the initial flux per correlation length is decreased, there is a point at which 
infinite strings disappear, as in the Hagedorn transition. 



I. INTRODUCTION 



Formation of linear defects (strings) at a symmetry-breaking phase transition is of great 
interest both in cosmology and in condensed matter physics. For a global symmetry break- 
ing, defects are formed because the phases of the Higgs field (order parameter) are uncor- 
related on scales greater than the characteristic correlation length ^h- In the simplest case 
of a global U{1) symmetry, a string is formed whenever thephase changes by 27r around a 
closed loop. This is the familiar Kibble-Zurek mechanism |2|. The statistical properties 
of the resulting string networks have been studied in numerical simulations [3|, with the 
following conclusions. 

The string network consists of two components: infinite strings and closed loops. The 
infinite strings constitute about 80% of the total string length and have the shape of random 
walks; their fractal dimension is Z) = 2, within statistical errors. The number of closed loops 
of length between / and / -|- dl per unit volume is given by 

n{l)dl = Ai~^'^l-^''^dl, (1) 

where A ~ 1 is a numerical coefficient. Long Brownian loops of length / have size R ~ 
(^i^/)^/^, and Eq. ([1]) gives a scale- invariant size distribution, 

n{R)dR ~ R-^dR. (2) 

Later work showed that these properties are very robust. In particular, simulations of Z2- 
string formation, which require a different symmetry breaking scheme, yield very similar 
results, except the fraction of length in closed loops is reduced from ~ 20% to ~ 6% 0, [sl- 
In the case of a gauge symmetry breaking, each string carries a quantum of magnetic 
gauge flux, and apart from the Kibble-Zurek mechanism, an additional mechanism of string 
formation can operate [6|. Any magnetic gauge field present before the phase transition will 
tend to be squeezed into quantized fiux tubes after the phase transition. This mechanism 
may operate in superconductors, where the stochastic magnetic field can be produced by 
thermal fiuctuations, and in cosmological phase transitions, where the field can be due either 
to thermal or to quantum fiuctuations. Flux trapping becomes the dominant mechanism 
of string formation when the magnetic field fiuctuations get sufficiently large, so that the 
typical area over which the magnetic fiux is equal to one fiux quantum is smaller than 
^"h 01 ■ This condition is often satisfied in superconductors and may well be satisfied in a 
cosmological setting. 

Another interesting example where fiux trapping may be important is the brane infiation 
model. Infiation in this model is driven by the attractive interaction energy of a and 
an anti-Z}3 brane separated in extra dimensions The branes eventually collide and 
annihilate, and cosmic Z)-strings (L)l-branes) can be produced in the process 



It has been argued in [101] that D-string formation by Kibble-Zurek mechanism is strongly 
suppressed in this model. However, a string network may still be formed by fiux trapping. 
Quantum fiuctuations of the gauge fields living on the branes get strongly amplified in the 
process of brane annihilation [l2[]. The resulting magnetic fields can then be squeezed into 
quantized fiux tubes to become D-strings. Fundamental {F) strings can be formed in a 
similar manner, by squeezing the electric component of another gauge field (orthogonal to 
that responsible for D-strings) into electric fiux tubes. ^ 



^ If both F and D strings are present, they can form bound states and combine into an FD network lCilll| 
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String networks formed by flux trapping may be rather different from those formed by 
the Kibble-Zurek mechanism. An important parameter here is the rms magnetic flux $c 
through an area where is the correlation length of the magnetic field. If $c is much 
greater than one flux quantum, then the strings will form in bundles containing many strings 
each [gI, 0, 14, 15]. This phenomenon has been observed in 2-dimensional simulations. 



where vortices (and anti-vortices) had a tendency to bundle together for large values of $c 
However, the properties of 3-dimensional string networks formed by flux trapping remain 
largely unknown. What are the shapes of long strings? Are they still Brownian? What 
fraction of the total length is in closed loops and what is their size distribution? If long 
strings form bundles, do they tend to stay in the same bundle, or switch from one bundle 
to another? To address these questions, we have developed a numerical simulation of string 
formation by flux trapping. 

Simulations of detailed Higgs and gauge field dynamics in a phase transition are compu- 
tationally expensive, and the maximum simulation size that we could achieve in this way 
would not be sufficient to get a good handle on the statistical properties of string networks. 
We therefore took a different approach. Our starting point is a stochastic magnetic field, 
which has presumably originated from either thermal or quantum fluctuations. We put the 
field on a lattice and use a relaxation technique to relax it to a state in which the flux 
through each plaquette is an integer multiple of the flux quantum. Strings can then be 
easily traced by following the flux lines from one cell to another. The net flux into each 
lattice cell remains equal to zero throughout the relaxation process, ensuring the continuity 
of strings. A "practical" application of these simulations is that they can be used as initial 
conditions for dynamical simulations of evolving string networks. 



II. SIMULATION 

A. Initial Conditions 

We have performed a series of numerical simulations borrowing techniques previously 
used in lattice field theory (see for example Ref. [l6|). We set up an initial configuration 
for the magnetic field on a cubic lattice by assigning the values of the vector field A(x) 
to the lattice links. It is then straightforward to compute the magnetic field flux through 
an individual plaquette from the values of the vector field on its boundaries. Since we are 
only interested in the magnetic part of the gauge field, it is sufficient to consider a vector 
potential of the form 

A(x) = 5^(2Vc.k)-^/'a(k)e^^■^ (3) 

k 

where Uk = \]<.\, V = is the volume of the simulation box, and we use the system of units 
in which h = c = 1. The reality of the electromagnetic field is ensured by the constraint, 
a(-k) = a*(k). 

We impose periodic boundary conditions, so the wave vector k takes a discrete set of 
values, k= (27r/L)n= (27r/L)(ni, ?t,2, 'n.3) with rii, ?t,2, ?t-3 = 0, ibl, .... The number of lattice 
points along each side of the box is = L/Ax, where Aa; is the lattice spacing. Wavelengths 
shorter than 2Ax cannot be represented in such a lattice, so we cut off the summation in 
(E]) at rij = ±N/2, 
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A(X^) = ^(47rLV|)-l/2jj(^)g2.*n.m/7V_ 
n 

Here, we have denoted the points on the lattice by = Ax m = Ax (mi, m2, m^). 

We obtain particular realizations of the vector field by drawing the values of the coef- 
ficients a(k) from an appropriate Gaussian distribution. In the case of superconductors, 
the distribution is expected to have a thermal (Rayleigh- Jeans) form, with a cutoff at short 
wavelengths [H], 

(a*(k)a,(k')) = -e-('=/'==)'%5kk', (5) 

where kc is related to the dissipation rate of the magnetic field, and K gives the amplitude 
of the spectrum, which in the thermal case is just the temperature. Defining the cutoff 
wavelength Ac = 2n/kc, we can rewrite this as 

k*(n)%(nO) = ^e"^'^'"'/'^'^^.^nn'. (6) 

This form of the spectrum for gauge field fiuctuations has also been found in brane 
infiation models, with fee ~ ~ I/t, where r is the characteristic timescale of brane 
annihilation [12i] . 



B. Relaxation technique 

Once the initial conditions for the magnetic field have been set, we need a mechanism 
that confines the field into strings carrying a fundamental unit of fiux, 

$0 = = a''/', (7) 

where e is the gauge coupling. Our final state should then have the magnetic fiux through 
all plaquettes equal either to zero or to an integral multiple of $o- This will allow us to 
follow the fiux around the simulation box to find the corresponding network of strings. 

In Appendix |A1 we compute the the rms magnetic fiux $c through a disk of diameter Ac 
in the spectrum of Eq. ([6]). We define the rms number of fiux quanta in such a disk, 

TV" = <l>c/<I>o ^ OAQ{KaX,y/^ . (8) 

For TV > 1, strings are expected to form bundles with several strings per bundle. The typical 
number of strings in a bundle is larger than Af, because bundles appear not in random places 
but where the fiux in a given direction is larger than usual. 

As mentioned before, the vector field A lives naturally on the links of the lattice and 
enters the computation of the magnetic fiux through the adjacent plaquettes. It is clear 
that changing the value of the gauge field on a link separating two plaquettes will only 
displace part of the fiux from one of those plaquettes to the other, keeping the total fiux 
unchanged. We would like to implement an algorithm that will drive the magnetic fiux 
through each plaquette to one of the values $ = n^o with n = 0, ±1, .... We accomplish this 
by evolving the vector field Ai on each link I according to the following equation of motion, 

dt ^ dAi ' ^ ' 
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where the sum is performed over all four plaquettes that have / as a side and V{^f) is a 
function of the magnetic flux through the f-th plaquette that has its minima at $ = n^Q. 
(It is assumed that V{—^) = V{^), so V{^) is independent of how one chooses the direction 
normal to the plaquette.) Eq. describes an overdamped system which is driven to the 
minimum of the potential energy, 



where the summation is taken over all the plaquettes in the lattice. The form of the function 
V{^) that we used in the simulation is given in Appendix [Bl This particular form is not 
essential for the method to work, although it is important to make sure that the final flux 
through a plaquette is approached in a smooth and monotonic way. This is achieved by 
choosing the time step At of the simulation sufficiently small. 

We finally have to choose a suitable lattice spacing Ax. A natural choice is to set it equal 
to the correlation length. Ax = A^;, as it was done in earlier simulations of string formation 
by the Kibble mechanism [1, 0, Hj. This works fine when the parameter A/" = $c/*^'o ~ 1, 
but for large values of J\f there is a large number of strings per lattice cell, which makes it 
hard to resolve individual strings. Hence, we used 



The overdamped evolution of Eq. converges on a state where the "potential energy" 
of Eq. ffTOl) vanishes because the fiux is properly quantized. However, this evolution may be 
quite slow, because some plaquettes may have fiux nearly equidistant between two multiples 
of $0 and so have little force driving them toward one or another. 

At a phase transition that confines fiux into strings, a loop of fiux can also shrink down 
to nothing. (For example, if there is a loop with less than half a quantum of fiux, this 
must happen.) In general, in a region with no net fiux to the outside (as we have because 
of periodic boundary conditions), in a slow phase transition the fiux will tend to diffuse to 
nothing, whereas in a rapid transition it will tend to be quantized. In the simulation, this 
choice is affected by the size of Ax. If Ax is taken so small that the average number of flux 
quanta through each plaquette is much less than 1, then the algorithm takes a long time to 
assemble a whole quantum of flux, and the flux tends to dissipate, so that there are very few 
strings. Thus even if Af is large, if we make Ax small we get primarily dissipation. Since 
we are interested in studying the shape of the string network, we will not make this choice. 
However, when A/" is small, even if we choose Ax ~ Ac as in Eq. (fTTl) . the average flux per 
plaquette will still be small, and so the diffusion of the flux may become important. 

III. RESULTS 
A. Large M 

We are mainly interested in studying the properties of a string network with a reasonably 
large number of strings per bundle. We have generated initial conditions for this type of 
conflguration by setting Ac = 6Ax and Af ~ 3.2. The result of the relaxation procedure 
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FIG. 1: Two-dimensional slice of a 64'^ lattice simulation. The black and white circles represent 
the position of vortices and anti- vortices respectively. Clusters of vortices (anti-vortices) are clearly 
visible in this figure, indicating that the corresponding strings in the 3D lattice form bundles that 
pierce this surface with the same orientation. We encircled the particular group of vortices that 
we study in detail in Fig. [2l 



shows that a typical 2D slice of the final configuration has several bunches, some of them 
including 10 strings or more. We show in Fig. [T]a 2D slice of a 64^ lattice simulation As one 



might expect, it resembles the 2D simulations discussed in Refs. |13l. |14J . 

To find out how long the string bunches stay together, we followed the three-dimensional 
trajectories of strings for a few steps in both directions, starting from the bunch encircled 
in Fig. [H The result is plotted in Fig. [2l We see that the strings only stay together for a 
few steps, quickly branching off in different directions (and possibly joining other bunches). 

We have also computed the fractal dimension of strings by averaging the distance R 
between points separated by length / along the string. We have computed this function R{1) 
by averaging over all sufficiently long strings present in 30 realizations; the result is shown 
in Fig. [31 

The best fit to the data is given by 

R = Af, (12) 

with 

A = 2.00 ±0.05 /? = 0.52 ±0.01, (13) 

where all lengths are measured in units of Ax. The stated errors here and below are an 
indication of the uncertainty in the parameter values, based on differences between separate 
runs. Equation (|T2|) is consistent with a random walk shape, /3 = 1/2, of the strings at large 
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FIG. 2: 3D structure of 5 of the 9 strings encircled in Fig. [TJ We also plot a small portion of the 
2D slice to indicate the region of space where the strings come together as they go through that 
plane. 
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FIG. 3: Average displacement between two points on the string versus the corresponding length 
along the string for the case with N ~ 3.2. Here and below white circles indicate points which are 
not included in the numerical fit. 

R. The effective step of the random walk, ^ ~ ~ 4, is comparable to Ac, as expected. It 
is also comparable to the characteristic length for which strings in a bundle stay together. 

Since our simulations were performed in a box with periodic boundary conditions, all 
the strings found in the lattice are, in fact, closed loops. The infinite string component of 
the network is represented by very long strings, which wind many times around the box. 
We counted all strings shorter than 2L as loops and the rest as infinite strings, but our 
results are insensitive to the particular choice of 2L as the separation between these two 
components of the network. We found that the fraction of energy stored in infinite strings 
is about 98% versus 2% in closed loops. 

The length distribution of loops is shown in Fig. |H The best fit to the data for sufficiently 
large loop sizes is 

n{l)dl = Bl-^'dl, (14) 

with the following parameters, 

5 = 0.017 ±0.001 7 = 2.50 ±0.03 (15) 

Once again, this is consistent with the scale- invariant distribution of Eqs. (ITPj) . 

We note that a string network with somewhat similar properties has been obtained by 
Pogosian and Vachaspati from a very different kind of simulation [l3|. They used a ran- 
dom phase distribution on a lattice, as in the Vachaspati- Vilenkin (VV) method but 
introduced additional phase differences 27in along the links, with the integer n taken from 
a Gaussian distribution. If the rms value of n is large, rirms 3> 1, then the typical number 
of strings through a plaquette is also large, and one can expect a 2D slice of the lattice 
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FIG. 4: Spectrum of closed string loops as a function of their length for the case with ~ 3.2. 

to exhibit bunches, with about rirms strings per plaquette. The number of closed loops is 
reduced compared to the standard VV simulation, because when a string returns to the 
same cell it started from, it does not necessarily close back on itself. Instead, it can now be 
connected to any one of ~ n^ms strings passing through the cell. 

B. M r^l 



To simulate a network with A/" ~ 1, we have chosen the following parameters: Ac = Ax, 
M ~ 0.28. These values have been selected to adjust the rms magnetic flux through a 
plaquette to 0.3$o- This roughly corresponds to the probabihty of ~ 30% of having a string 
through a plaquette in Kibble mechanism simulations of With this set of parameters, 
the magnetic fluxes through neighboring plaquettes are only weakly correlated with one 
another, so we expect that strings will not form bundles and that the resulting network 
will be similar to the ones obtained in And indeed, the results we obtained are nearly 
identical to those of , even though the numerical methods used in the two simulations are 
rather different in nature. 

We have performed these simulations in cubic lattices of size 16'^, 32'^ and 64^ and found 
that the average fraction of string length in long (infinite) strings was respectively 86%, 84% 
and 82%. These values are very close to those found in Ref. |3|. 

The distance versus length along the string, -R(/), and the length distribution of loops, 
n(/), are plotted in Figs. and [HI respectively, for that largest box of size 64^. The data are 
well fitted by Eqs. ([12]) , ([H]) with 



A = 1.45 ±0.03 /3 = 0.51 ±0.01, 



(16) 
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FIG. 5: Average displacement between two points on the string versus the corresponding length 
along the string for the case with M ~ 0.28. 
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FIG. 6: Spectrum of closed string loops as a function of their length for the case with A/" ~ 0.28. 
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FIG. 7: Fraction of the total string length in infinite strings as a function of the parameter M . 



B = 0.065 ±0.003 7 = 2.51 ±0.02. 



(17) 



As before, the strings have Brownian shapes and the loops exhibit a scale-invariant distri- 
bution. 



C. Small M 

To explore the properties of the network for A/" ^ 1, we kept Ac = Ax and gradually 
decreased A/", starting from the value of 0.28 cited in the preceding subsection. For A/" <^ 1, 
the rms flux through a plaquette is much smaller than $o- As a result, most plaquettes have 
no strings. The strings have the form of small closed loops and appear in rare places where 
the magnetic field fluctuates well above its rms value. 

In Fig. [7] we plot the fraction of length in infinite strings as a function of M . The 
simulations discussed in subsections III. A and III.B above correspond to A/" = 3.2 and 
M = 0.28, respectively. As A^ is decreased, infinite strings constitute a smaller and smaller 
fraction of the total string length and finally disappear at AC: ~ 0.15. The total string 
energy density also decreases with the decrease of A^ (see Fig. [8]), but remains finite at (and 
below) the critical value A^. For all values of Af > Afc the loop distribution retains its 
scale- invariant form ([1]). For A/" < A^, most of the string length is in the smallest loops, and 
the number density of loops n{l) decreases rapidly with the loop's length /. 

It is interesting to note that very similar behavior is observed in a thermal ensemble of 
strings with a lower cutoff on the loop length, Imin- As the string length per unit volume 
is decreased below a certain critical value, pc ~ the system undergoes the Hagedorn 
transition, characterized by the disappearance of infinite strings. It was noted earlier [18, 1^ 
that string configurations resulting from simulations of Kibble-Zurek-type phase transitions 
resemble those obtained from a thermal string ensemble. 
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FIG. 8: Plot of the average total string length (black circles) as well as the different contribu- 
tions from infinite (squares) and closed strings (white circles) as a function of M. The average is 
performed over 100 simulations of a 32^ lattice. 

However, we also note that in this regime, the average flux per plaquette is quite small. 
This appears to increase the effects of diffusion of the flux as compared to confinement, and 
so it is possible that the decrease in long strings is due to the increased diffusion. 

IV. CONCLUSIONS 



We have studied the statistical properties of string networks formed by the flux trapping 
mechanism. As anticipated in earlier analyses 0, [l3, 14, 15j], the character of the network de- 
pends on the rms flux through a correlation- size region, $c, or equivalently on the parameter 
Af = ^c/^o, where $o is the flux quantum. 

We found that for large values of A/" most of the string length is in a tangled network 
of infinite strings, having the shape of random walks, with closed loops contributing only 
a small fraction of the total length. The step of the random walk is comparable to the 
correlation length Ac of the magnetic field. The typical inter-string separation, d ~ Ac/ 
can be much smaller. A 2D slice of the simulation exhibits bundles of strings pointing in the 
same direction, with ~ A/" strings per bundle. We found, however, that strings belonging 
to the same bundle do not stay together for long. As we follow them in 3D, they part each 
other's company on a length scale not much exceeding Ac. 

For A/" ~ 1, there are no string bundles, and the properties of the string network are 
nearly identical to those of the networks formed by the Kibble-Zurek mechanism. 

For small values of Af, the string density is exponentially suppressed and most of the 
string length is in the form of small loops of length / ~ Ac. 

How does the initial string configuration affect the subsequent evolution of the network? 
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For TV ^ 1, the properties of the network differ from the standard initial configuration of 
Ref. |3i] on relatively short length scales, N'~^^'^\c ^ I \c, which can have an effect on 
early string evolution. These effects can be studied by dynamical string simulations in an 
expanding universe, with initial conditions generated by simulations of the type described 
here. At late times, the string bundles break apart, and we expect the network evolution 



to follow the standard scenario (see, e.g., |20|]). For type-I strings, which can have arbitrary 



large winding numbers, the initial configuration with Af ^ 1 may facilitate the formation of 
higher- winding strings^. This may have an effect on the subsequent evolution of the network. 
For A/" ~ 1, the initial configuration is similar to the standard one, so the standard evolution 
scenario is followed from the very beginning. Finally, for A/" ^ 1, all strings are in the form 
of small closed loops, which rapidly shrink and disintegrate. 
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APPENDIX A: AVERAGE FLUX THROUGH A DISK 

We want to compute the rms magnetic flux through a disc of radius R. We follow similar 



calculations done in [IJ, Il5[. The mean square flux is 

= (^j\^xd^yB,ix,z = 0)B,iy,z = 0)^ . (Al) 
Using our definition for A given in Eq. ([3]), this becomes 

mnf) = r d^xd^yS^iW^u^u^)^-^/^^ [g,fc,.(a;(q)a,(k)) + g,fc,(a:.(q)a,(k))] e^(k— 'i-y) . 

By X and y we mean the 3-vectors whose z component is 0, and we have neglected the cross 
terms that average to zero. Using now equation (jS]) we arrive at 

($(i?)2) = r d''xd^yy{2Vu^)-\kl + P)— e-(^/^=)'e^'^-('^-y) , (A2) 

where k is the length of k. We now take the infinite volume limit and replace the sum by 
an integral, using 

^E^T^/^'^' (A3) 



so we can rewrite the expression above in the following form. 



($(i?)2) = _1_J^ d'xd'y J d3fc|±.e-('^/'==)'e*^ (^-y) , (A4) 



^ Field theory simulations ^] indicate that this is indeed the case, although much bigger simulations are 
needed in order to obtain string bunches with A/" ^ 1. 
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where we have defined k± — {k^ + kyY/"^. We can exchange the order of the integrals to get, 



Performing the integral in r we obtain 

KR^ fd^k 



mR)') = ^e-^'/'^^^' Mk^Rf (A5) 

We now break up the integral into an integral in k^, a radial integral in k±, and a trivial 
angular integral, to get 

($(i?)2) = y dk^k^Ji{k±R)^e-^'/'^^ J dk 



?2 r°o 



dk±Ji{k±R)^eTic{k/kc) = ^i— ^ 



U2 I U2 

1 



where 



^W = /„ 5.5.3.--) (A6) 

where is a hypergeometric function. 

We can now define J\f as the rms fiux through a disk of diameter Ac, i.e., R = Ac/2, in 
fundamental fiux units, namely, 

^ = (^(V2)T/^ ^ (X«Ac) V2 (^^(2)) 0A6{KaX^y/' (A7) 



APPENDIX B: THE CONFINING POTENTIAL V{<^>). 

We wish to choose a potential which will confine the field by driving the fiux through 

each plaquette to an integer multiple of $o- We treat any multiple the same as any other, 
so V should be periodic with period $o- Similarly, upward and downward fiux are to be 
treated the same, so V{^) = V{-^). 

We choose an inverted parabola centered at $o/2, so that the fiux will be driven to the 
nearest multiple of $o- If we continue this parabola all the way to zero, we find that the 
fiux jumps back and forth around zero, instead of settling smoothly to zero. Thus we use 
a smooth quadratic form for near zero. The smoothing extends to a critical value 

= Any larger value can result in one unit fiux being spread out throughout the 

lattice instead of being concentrated on a single plaquette. 

Our potential is then given by 

vm = H . (Bi) 

^ ^ I 1 - C2 ($0/2 - < $ < $0/2 ^ ^ 
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FIG. 9: Confining potential as a function of the magnetic flux through a single plaquette. 



and for other values of $ by reflection and periodicity, with 

2 



Cl 



2 



C2 



(B2) 



(B3) 



<l>o($o/2-$*) ' 

which make V{^) and its derivative continuous at $*. We plot in Fig. [9] the potential V{^) 
for N = 64. 

We then choose the largest time step At which will not allow the potential to overshoot 
in a single step. 
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